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I report on recent work done in collaboration with Irinel Caprini and Gilberto Colangelo 0. We 
observe that the Roy equations lead to a representation of the tt-k scattering amplitude that ex- 
clusively involves observable quantities, but is valid for complex values of s. At low energies, this 
representation is dominated by the contributions from the two subtraction constants, which are 
known to remarkable precision from the low energy theorems of chiral perturbation theory. Evalu- 
ating the remaining contributions on the basis of the available data, we demonstrate that the lowest 
resonance carries the quantum numbers of the vacuum and occurs in the vicinity of the threshold. 
Although the uncertainties in the data are substantial, the pole position can be calculated quite 
accurately, because it occurs in the region where the amplitude is dominated by the subtractions. 
The calculation neatly illustrates the fact that the dynamics of the Goldstone bosons is governed 
by the symmetries of QCD. 

PACS numbers: 11.30.Rd, 11.55.Fv, 11.80.Et, 12.39.Fe, 13.75.Lb 



Pions play a crucial role whenever the strong interac- 
tion is involved at low energies - the Standard Model pre- 
diction for the muon magnetic moment provides a good 
illustration. The present talk concerns the remarkable 
theoretical progress made in low energy pion physics in 
recent years. I concentrate on work based on disper- 
sion theory Q, but also show some lattice results, ob- 
tained with quark masses that are small enough for a 
controlled extrapolation to the values of physical interest 
to be within reach. The low energy precision measure- 
ments 0, 0, Q| are consistent with the predictions and 
even offer a stringent test for one of these. 

From the point of view of dispersion theory, 7T7t scat- 
tering is particularly simple: the s-, t- and u-channels 
represent the same physical process. As a consequence, 
the real part of the scattering amplitude can be repre- 
sented as a dispersion integral over the imaginary part 
and the integral exclusively extends over the physical 
region 5j. The representation involves two subtraction 
constants, which may be identified with the S-wave scat- 
tering lengths Oq, a^. The projection of the amplitude on 
the partial waves leads to a dispersive representation for 
these, the Roy equations. 

The pioneering work on the physics of the Roy equa- 
tions was carried out more than 30 years ago [fj]. The 
main problem encountered at that time was that the two 
subtraction constants occurring in these equations were 
not known. These constants dominate the dispersive rep- 
resentation at low energies, but since the data available 
at the time were consistent with a very broad range of 
S- wave scattering lengths, the Roy equation analysis was 
not conclusive. 

The insights gained by means of xPT thoroughly 
changed the situation. The corrections to Weinberg's 
low energy theorems pj for af], ag (left dot in Fig.l) have 
been worked out to first non-leading order [8| (middle 



dot) and those of next-to-next-to leading order are also 
known 9] (dot on the right). Very accurate predictions 
for the scattering lengths are obtained by matching the 
chiral and dispersive representations in the interior of 
the Mandelstam triangle [lCj (small ellipse). The lattice 
results of the MILC collaboration also yield an esti- 
mate for the scattering lengths. Using their values for the 
coupling constants L4, L5, L§, Lg of the effective chiral la- 
grangian and neglecting two loop effects, we arrive at the 
one standard deviation contour indicated by the small 
dashed ellipse. The horizontal band represents the value 
of the scattering length a§ obtained by the NPLQCD 
collaboration 

The plot shows that 7T7t scattering is one of the very 
rare hadronic processes where theory is ahead of experi- 
ment: the two large ellipses and the tilted band indicate 
the results obtained on the basis of experiments done at 
Brookhaven and CERN 0, 3 ■ Possibly, the error es- 
timates attached to the lattice results are too optimistic 
and the analysis of the K — > 3tt data also requires further 
study, but it is fair to say that all of the published results 
confirm the theoretical predictions. 

In the following, I focus on the results obtained on 
this basis for the low energy properties of the isoscalar 
S'-wave. The corresponding S'-matrix element, S® = 
77$ exp 2 z <5g, is related to the partial wave amplitude 
by 



S$(s) = l + 2ip( S )t°(s), p(s) = Vl-4M2/s. (1) 

The various phenomenological analyses are not in good 
agreement 13]. In fact, until ten years ago, the infor- 
mation about the a — Jo (600) was so shaky that this 
resonance was banned from the data tables. The work of 
Tornqvist and Roos |l4( resurrected it, but the estimate 
M CT -ir CT = (400-=- 1200) -i (300-^500) MeV of the Par- 
ticle Data Group |l5j | indicates that it is not even known 
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FIG. 1: 5* wave scattering lengths. 



for sure whether the lowest resonance of QCD carries the 
quantum numbers of the a or those of the p. 

The positions of the poles represent universal prop- 
erties of the strong interaction which are unambiguous 
even if the width of the resonance turns out to be large, 
but they concern the non-perturbative domain, where an 
analysis in terms of the local degrees of freedom of QCD 
- quarks and gluons - is not in sight. One of the rea- 
sons why the values for the pole position of the a quoted 
by the Particle Data Group cover a very broad range is 
that all of these rely on the extrapolation of hand made 
parametrizations: the data are represented in terms of 
suitable functions on the real axis and the position of the 
pole is determined by continuing this representation into 
the complex plane. If the width of the resonance is small, 
the ambiguities inherent in the choice of the parametriza- 
tion do not significantly affect the result, but the width 
of the a is not small. 

We have found a method that does not invoke 
parametrizations of the data. It relies on the following 
two observations: 



1. The S'-matrix has a pole on the second sheet if 
and only if it has a zero on the first sheet. In order 
to determine the pole position it thus suffices to have a 
reliable representation of the scattering amplitude on the 
first sheet. 

2. The Roy equations hold not only on the real axis, 
but in a limited region of the first sheet. Since the pole 
from the a occurs in that region, we do not need to invent 
a parametrization, but can rely on the explicit represen- 
tation of the amplitude provided by these equations. 

Roy established the validity of his equations from first 
principles for real values of s in the interval — 4M^ < s < 
60M^ . Using known results of general quantum field the- 
ory |l6( , we have demonstrated that these equations also 
hold for complex values of s, in the intersection of the 
relevant Lehmann-Martin ellipses lj. The dash-dotted 
curve in Fig |3 shows the domain of validity that follows 
from axiomatic field theory, while the full line depicts 
the slightly larger domain obtained under the assump- 
tion that the scattering amplitude obeys the Mandelstam 
representation. I emphasize that the boundary does not 



3 



40 




-40 1 1 1 1 ' 

20 40 60 

2 

Res in units of M 

n 



FIG. 2: Domain of validity of the Roy equations. 



represent a singularity of the amplitude, but merely lim- 
its the region where the Roy equations hold in the form 
given. Modified representations with a much larger do- 
main of validity can be found in the work of Roy and 
Wanders .17] and in the references quoted therein. 

For our analysis, it is essential that the dispersion in- 
tegrals are dominated by the contributions from the low 
energy region: because the Roy equations involve two 
subtractions, the kernels fall off with the third power 
of the variable of integration. The left hand cut plays 
an important role here: taken by itself, the contribution 
from the right hand cut is sensitive to the poorly known 
high energy behaviour of the imaginary parts, but taken 
together with the one from the left hand cut, the high 
energy tails cancel. In this connection, I note that most 
pole determinations assume that the left hand cut can be 
neglected |l8|. Since the distance between the pole and 
the left hand cut is not large, this assumption is difficult 
to justify. 

The Roy equations thus provide us with an explicit 
representation of the function Sq(s) for complex values 
of s, in terms of rapidly convergent dispersion integrals 
over the imaginary parts of the partial waves. In con- 
nection with the determination of the pole from the a, 
the most important contribution is the one from the sub- 
traction term. The dispersion integrals over the S- P- 
and D-waves generate a correction which can be evalu- 
ated with available phase shift analyses - in particular 
with the one obtained by solving the Roy equations ^(| . 
The contributions from high energies and high angular 
momenta can be estimated by means of the Regge rep- 
resentation of the scattering amplitude [2(j - these con- 
tributions barely affect the pole position. 

For the central solution of the Roy equations, the func- 
tion Sq(s) contains two pairs of zeros in the domain of 
interest: 

s a = (6.2 ± U2.3) Ml , s /o = (51.4 ±i 1.4) M^. (2) 



These are indicated in Fig0 which may also be viewed 
as a picture of the second sheet - the dots then repre- 
sent poles rather than zeros. For comparison, the figure 
also indicates the position of the zeros in Sl(s), which 
characterize the p. 

The higher one of the two pairs of zeros represents the 
well-established resonance Jo (980), which sits close to the 
threshold of the transition 7T7t — > KK. The correspond- 
ing pole generates a spectacular interference phenomenon 
with the branch point singularity, which gives rise to a 
sharp drop in the elasticity. Our analysis adds little to 
the detailed knowledge of that structure. 

The lower pair of zeros corresponds to a pole in the 
lower half of the second sheet at [lj 

M a - l - T a = ^ = 441 - i 272 + J 12 . 5 MeV . (3) 

The error bars account for all sources of uncertainty. 
They are calculated by (a) estimating the uncertainties 
in the input used when solving the Roy equations and (b) 
following error propagation to determine the uncertainty 
in the result for the pole position. For details, I refer to 
our publication |l|. 
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FIG. 3: Dominance of the subtraction term. 

The reason why the pole position of the a can be cal- 
culated rather accurately is that (a) the pole occurs at 
low energies and (b) there, the isoscalar 5*-wave is dom- 
inated by the subtraction term. This is illustrated in 
FigEl where the real part of t%(s) is plotted versus s. 
The graph demonstrates that, in the region shown, the 
full amplitude closely follows the subtraction term. The 
final state interaction does generate curvature - in par- 
ticular, the cusp at s = AMI i s visible - but since Gold- 
stone bosons of low momentum interact only weakly, the 
contributions from the dispersion integrals amount to a 
small correction. In the language of xPT , the disper- 
sion integrals only show up at NLO. Moreover, at lead- 
ing order, the subtraction constants are determined by 



4 



the pion decay constant. Dropping the dispersion inte- 
grals and inserting the lowest order predictions for the 
scattering lengths, the Roy equation for the isoscalar S- 
wave reduces to the well-known formula, which Weinberg 
derived 40 years ago @, 



t° (s) 



2s- Ml 
32ttF2 



(4) 



and which is shown as a dash-dotted line in Fig|3| 

The main feature at low energies is the occurrence of an 
Adler zero, t%{sA) — 0. At LO, the zero is at sa = \M%. 
The higher order corrections generate a small shift, which 
can be evaluated from our solution of the Roy equations. 
The uncertainties in the result, sa — (0.41±0.06) M%, are 
dominated by those in our predictions for the scattering 
lengths. While the Adler zero sits below the threshold, 
at a real value of s, the lowest zeros of the S'-matrix, 
So(s a ) = 0, occur in the region Res a > 4M 2 , with an 
imaginary part that is about twice as large as the real 
part. 

The formula |J3J explains why the S'-matrix has a zero 
in the vicinity of the threshold. In this approximation, 
the zero of 5g(s) occurs at = 365 - i 291 MeV. The 
number differs from the "exact" result in Eq. by about 
20 percent. The dispersion integrals are essential for the 
partial wave to obey unitarity, but only represent a cor- 
rection. In view of this, the precision in Eq. J3J is rather 
modest. 

I conclude that the same theoretical framework that 
leads to incredibly sharp predictions for the threshold 
parameters of tttt scattering [10J also shows that the low- 
est resonance of QCD carries the quantum numbers of 
the vacuum. The physics of the a is governed by the 
dynamics of the Goldstone bosons: the properties of the 
interaction among two pions are relevant 21] . In quark 
model language, the wave function contains an important 
tetra-quark component p^ . 

The resonance /o(980) is also governed by Goldstone 
boson dynamics - two kaons in that case. Very recently, 
the method described above was applied to the case of 
irK scattering [2*3 . In this case, the analogue of the 
back-of-the-envelope calculation sketched above relies on 
the tree level approximation for the I = h S-w&ve ob- 
tained from the effective SU(3)xSU(3) lagrangian and 
yields m K = 671 — i 262 MeV, remarkably close to the 
"exact" value, obtained from the solution of the Roy- 
Steiner equations, m K = (658 ±13) -i (278.5 ± 12) MeV. 
Evidently, the physics of the k is very similar to the one 
of the a. 

The BES data [24| on the decay J/ip — > lotttt play a 
prominent role in scalar meson spectroscopy. Although I 
did not discuss this at the conference, I add a comment 
concerning the recent claim [25j that our analysis is in 
conflict with these data. I denote the projection of the 
decay amplitude onto the tttt configuration with I — £ — 
by Aq(s). If rescattering on the u> and the inelasticity 
due to 4-7T final states is neglected, unitarity implies that 
Aq(s) shares an important property with the scalar pion 



form factor Fq(s) = (n + ir~ out|uu |0): below the KK 
threshold, the phase coincides with the phase shift 6q(s) 
of -kit scattering j2f|. On the interval 4M 2 < s < 4M^, 
the ratio R(s) = Aq(s) / Fq(s) must therefore be real and 
can vary only slowly. 

FigEt shows the magnitude of this ratio, obtaine d by 
dividing the outcome of the partial wave analysis |27j 
for |An (s)| with the dispersive result for the scalar form 
factor |23 (arbitrary units on the vertical axis). The 
ratio indeed varies only slowly with the energy 29] . This 
demonstrates that the profile of the BES data on |v!q (s) | 
closely resembles the one of the scalar form factor. In 
either case, there is a peak around y/s ~ M a , but the 
phase reaches 90° only when the energy becomes about 
twice as large as the mass of the a. 

The low energy behaviour of the various amplitudes 
can be understood on the basis of twice subtracted dis- 
persion relations. In the case of nn scattering, the sub- 
traction constants are predicted by chiral symmetry. As 
a consequence, the low energy behaviour of the phase 
shift Sq(s) and the position of the pole from the a are 
understood. The symmetry, in particular, requires the 
occurrence of an Adler zero, which suppresses the scat- 
tering amplitude at low energies, so that the peak in that 
amplitude is pushed up, to yfs ~ 2 M a . There is no such 
suppression in the decay amplitude, nor is there a pre- 
diction in that case: the two subtraction constants in the 
dispersion relation for R(s) = Ao(s)/Fq(s), 



R(s) = Ro + RiS- 



( S -2Mi) 2 



dx Imi?(x) 



{x-2M 2 K f (x-s) 



must be taken from experiment. The integral receives a 
contribution from s > 4M|-, as well as one from s < 0. 
For definiteness, I have identified the subtraction point 
with the center of the interval between these two cuts. 
The straight line shown in Fig^Ji is obtained by setting 
so = —Rq/Ri = 1.65 GeV 2 and simply dropping the dis- 
persion integral. Fig^Ja shows that the approximation 
A (s) ~ Ri(s — s )F (s) indeed accounts for the pro- 
nounced energy dependence seen in the BES data. 

Above the range covered by the data points, the disper- 
sive contributions presumably become quite important. 
In particular, the behaviour near the KK threshold need 
not be the same as the one of the non-strange form fac- 
tor Fq(s). As can be seen in Fig0Ji, this form factor has 
a narrow dip there. The form factor of the operator ss 
instead exhibits a narrow peak 28]. If all final states 
other than 7T7t and KK are ignored, the general solution 
of the unitarity conditions is a linear combination of the 
two. An experimental investigation of the structure of 
A (s) in the region of the /o(980) could shed some light 
on the importance of strange quarks in this context. At 
any rate, the zero, which a linear approximation for the 
ratio R(s) necessarily entails, sits far beyond the region 
where that approximation is meaningful. 

I do not see any reason to doubt that the BES data 
and the partial wave analysis thereof are correct. These 



0.2 0.4 0.6 0.8 1 1.2 

GeV 2 



0.2 0.4 0.6 0.8 1 

GeV 



FIG. 4: Jip — > conn decay: energy dependence of the amplitude |^4o(s)|- 



clearly reveal the presence of the a, but the pole po- 
sitions extracted from simple Breit-Wigner treatments 
1^1 of the amplitude Aq(s) or more elaborate models 
|25| are subject to large theoretical uncertainties: the 
extrapolation off the real axis is sensitive to details of 
the parametrization that are not understood. In the lan- 
guage of dispersion theory, the problems arise from (a) 
the inelastic channels mr — > KK and tttt — > rjr], (b) the 
/o(980), (c) the left hand cut and (d) rescattering on the 
lu. Unitarity and crossing symmetry very strongly con- 
strain the scattering amplitude, but barely tell us any- 



thing about the transition J '/ip — » lout:. 

We have applied our method to the models for the scat- 
tering amplitude proposed by Bugg |25|. The outcome 
for the mass and width of the a is close to the central 
values in Eq. I conclude that our pole position is 
consistent with the BES data. There is a conflict only 
with the values for the mass and width obtained from the 
extrapolation of those models to complex values of s. It 
arises because the quoted errors do not account for the 
uncertainties inherent in such extrapolations |32| . 
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